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ABSTRACT 


With  the  availability  of  an  ELF  transmitter  (the  Bravo  test  facility 
at  Clam  Lake,  Wisconsin)  it  now  becomes  feasible  to  make  far-field  field 
strength  measurements  at  the  lower  ELF  frequencies.  Field  strength 
measurements  were  made  in  Hawaii  and  California  at  45  and  7  5  Hz  in  an 
attempt  to  determine  the  east-west  attenuation  rate  under  daytime  or  night¬ 
time  conditions  over  the  propagation  path.  As  a  result  of  high  atmospheric 
noise  conditions  and  a  low  transmitted  power  level,  the  relatively  large 
uncertainty  in  signal  estimates  makes  accurate  estimation  of  <y  impossible. 
However,  the  data  indicate  a  daytime  east  to  west  attenuation  rate 
between  0.7  and  1.2  db/Mm  and  a  nighttime  rate  between  1.0  and  1.8  db/Mm. 

Future  experiments  are  being  planned  to  refine  the  estimates  of  <y  to 
within  ±0.2  db/Mm. 


Accepted  for  the  Air  Force 

Joseph  R.  Waterman,  Lt.  Col.  ,  USAF 

Chief,  Lincoln  Laboratory  Project  Office 
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I. 


INTRODUCTION 


Over  a  three  week  period  covering  July  6  through  July  24  (1970), 
the  test  facility  transmitter  in  Wisconsin  was  tasked  to  transmit  a  sinusoidal 
signal  for  four  hour  intervals  during  which  the  Wisconsin  -  Hawaii  propaga¬ 
tion  path  was  either  entirely  in  daylight  or  in  nighttime  conditions.  Trans¬ 
missions  were  scheduled  from  either  0100  to  0500  or  1300  to  1700  hours 
Central  Daylight  Savings  Time.  For  the  first  two  weeks  transmissions  were 
scheduled  at  45Hz  and  for  the  third  week  at  75  Hz.  The  antenna  current 
level  was  150A*  .  Receivers  were  positioned  at  Hawaii's  Volcano  National 
Park  and  at  Hat  Creek,  California  along  a  great  circle  path  which  includes 
the  Wisconsin  transmitter  facility. 

Until  recently  theoretical  estimates  and  sferic  measurements  of  the 
attenuation  coefficient  (ck)  at  45  Hz  resulted  in  a  range  of  values  between  0.6 
and  0.  8  db/Mm  which  were  assumed  applicable  under  nighttime  as  well  as 
daytime  conditions  and  were  virtually  independent  of  propagation  direction. 
Recent  theoretical  considerations,  which  incorporate  the  effects  of  ions, 
have  suggested  nighttime  attenuation  rates  of  the  order  of  1.  1  to  1.2  db/Mm 
[  1,2]  .  These  recent  theoretical  results  still  indicate  only  a  very  small 
directional  dependence  (i.  e.  ,  ~  0.  1  db/Mm  difference  in  nighttime  attenuation 
coefficient  for  propagation  from  east  to  west  versus  that  of  west  to  east). 
However,  two  recent  station  sferic  measurements  by  Smith  [3]  and  Hughes  [4] 
have  been  interpreted  to  indicate  east  to  west  attenuation  coefficients  for 
daytime  conditions  which  are  greater  than  2. 0  db/Mm.  As  a  result  of  this 
present  uncertainty  in  knowledge  of  the  attenuation  coefficient's  diurnal  as 
well  as  directional  dependence,  it  becomes  difficult  to  decide  on  an 
appropriate  sizing  for  an  ELF  communications  system.  Clearly  a  series  of 
measurements  to  determine  the  attenuation  coefficient  within  a  few  tenths  of 
a  db/Mm  in  various  propagation  directions  and  for  different  ionospheric 
conditions  would  be  invaluable  in  making  an  economical  system  sizing. 

Except  for  the  nighttime  transmissions  on  July  9-10  and  July  10-11. 
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However,  for  the  present  radiated  power  levels  this  would  be  an  extremely 
difficult  measurement  for  reasons  which  will  be  explained  more  fully  in  the 
remainder  of  the  report.  With  the  experimental  configuration  as  set  up  in 
Hawaii  and  California  it  is  possible  to  determine  whether  an  east  to  west 
daytime  attenuation  coefficient  greater  than  2.  0  db/Mm  did  exist  as  claimed 
by  Hughes.  Accordingly  the  main  purpose  of  the  recently  completed  measure¬ 
ments  was  to  show  that  extremely  large  values  of  attenuation  coefficient  (in 
excess  of  2.0  db/Mm)  were  not  in  evidence.  A  subsidiary  purpose  was  to 
test  experimental  techniques  which  are  to  be  used  in  a  future  expanded  series 
of  measurements  which  will  attempt  to  measure  a  within  a  few  tenths  of  a 
db/Mm.  These  future  experiments  are  described  in  detail  in  another  report. 

It  should  be  recognized  that  a  is  just  one  of  the  parameters  in  the 
analytic  expression  for  field  strength  at  any  point  in  the  earth  ionosphere 
cavity.  The  analytic  expression  most  commonly  used  is  an  asymptotic  form 
of  the  exact  solution  of  propagation  in  a  cavity  bounded  by  two  concentric 
spherical  surfaces  representing  the  earth  and  the  "effective  fl  lower  boundary 
of  the  ionosphere.  Each  of  these  surfaces  is  assumed  to  have  a  uniform 
conductivity  over  its  entire  surface.  Moreover  the  transmitter  is  assumed 
to  be  a  short  vertical  dipole.  In  reality  however  the  earth Ts  conductivity 
varies  from  point  to  point  on  its  surface,  the  antenna  is  in  fact  a  horizontal 
line  and  the  ionosphere  changes  significantly  from  nighttime  to  daylight 
conditions  and  has  a  continuous  profile  with  height  rather  than  a  step  dis¬ 
continuity.  Rather  than  attempt  a  full  wave  solution  on  a  realistic  model  it 
has  been  relatively  fruitful  to  make  intuitive  modifications  to  the  ideal 
homogeneous  concentric  cavity  solution.  The  main  problem  of  course  is  to 
determine  the  magnitude  of  error  that  these  approximations  introduce.  Wait[5] 
provides  a  good  review  of  the  various  approximations  for  near  field,  far 
fields,  field  near  antipode  etc.  In  the  far  field  and  away  from  the  antipode, 
the  horizontal  magnetic  intensity  for  the  direct  component  traveling  along  the 
shorter  segment  of  the  great  circle  joining  the  transmitter  and  receiver  is 
given  by 
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angle  between  axis  of  horizontal  antenna  and  propagation  path 
direction 

ground  conductivity  at  transmitter 

effective  height  of  ionosphere 

phase  velocity  of  the  only  propagating  mode 
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earth  radius 
u>  /  c  =  2  77/  X 

distance  between  receiver  and  transmitter 
attenuation  coefficient 


Now  there  is  also  another  magnetic  field  component  which  travels  around  the 

world  along  the  longer  segment  of  the  same  great  circle  path.  This  component 

at  the  receiver  differs  from  equation  (1)  in  that  one  only  needs  replace  e  °  ^ 

-0? '(2ff r  -d)  jk  (c/V  ,)d  •>ko(c/VDh)(2lTre'd) 

by  e  and  e  ^  by  e  ”  .  (In  addition  a 

77/2  phase  shift  is  introduced  as  this  wave  passes  through  the  antipode.  )  The 

direct  and  the  around-the -world  components  add  to  give  a  standing  wave  pattern. 

The  positions  of  the  standing  wave  minima  and  maxima  depend  on  the  phases 

of  the  direct  and  round-the-world  wave  components.  For  a  realistic  global 

ionospheric  model,  with  local  (and  perhaps  directional)  changes  in  a  and 

V  ^  one  can  make  the  approximation 
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(D  refers  to  direct  wave  and  R  the  round-the-world  component.  ) 
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The  total  field  standing  wave  pattern  then  is  the  sum  of  the  above  two 

expressions.  The  main  point  is  that  if  one  desires  to  measure  the  parameter 

say  by  a  two  station  measurement,  it  is  imperative  to  know  the  phases 

of  both  (H  )r>  and  (H  )~  at  both  stations  in  order  to  determine  receiver 
(0  -K  <p  U 

location  with  respect  to  the  maxima  and  minima  of  the  standing  wave  pattern. 

With  the  phases  known,  one  still  must  have  a  reasonable  estimate  of 
Zirr  -d 

Of^dr  so  that  can  be  estimated  as  well.  It  should  be  evident 
"  o  R  D 

that  no  single  frequency,  two  fixed  receiver  measurements  can  fix  the  value 
of  Qfj-j  with  any  high  degree  of  accuracy  unless  the  other  parameters  are  known 
as  well.  However,  even  if  we  have  total  ignorance  of  the  values  of 
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The  uncertainty  in  the  structure  of  the  standing  wave  pattern,  however, 
is  only  one  element  that  produces  an  uncertainty  band  about  The  fact  is 

that  the  relatively  low  S/N  ratios  encountered  are  responsible  for  estimates 


4 


of  field  strengths  (and  hence  estimates  of  o^)  that  have  reasonably  large 
confidence  intervals  associated  with  the  signal  level  estimates.  These  two 
factors  then  are  the  prime  reasons  why  with  this  limited  measurement  series, 
we  can  do  little  more  than  determine  whether  a ^  is  greater  or  less  than 
2,0  db  /Mm. 

Part  II  will  briefly  describe  first  the  method  of  recording  data  and 
second  the  method  by  which  the  data  is  then  processed  and  finally  analyzed. 

Part  III  records  the  data  successfully  recorded  and  analyzed  and 
catalogs  the  results.  In  this  section  we  will  show  that  full  potentiality  of 
this  experiment  was  degraded  due  to  the  high  values  of  atmospheric  noise 
(and  as  yet  occasional  unidentified  interference  in  Hawaii).  The  experiment 
was  originally  planned  for  April  when  the  atmospheric  noise  levels  are  low 
(about  -135  dbH/ J  Hz  ).  Due  to  delays  in  transmitter  mitigation  the  higher 
current  levels  were  not  available  until  July.  In  July  the  late  morning  summer 
noise  levels  typically  were  5  to  10  db  higher  than  -135  dbH/^/Hz.  Hence  the 

A 

S/N  ratios  were  reduced  from  5  to  10  db  from  what  we  had  originally  expected. 
This  large  loss  in  S /N  significantly  increases  the  confidence  interval 
(uncertainty)  associated  with  the  signal  estimate. 

Part  IV  interprets  the  results  with  regard  to  answering  the  question 
of  how  big  the  east-to-west  attenuation  coefficient  is.  Part  V  is  a  short 
summary  of  the  principal  results. 
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II. 


DESCRIPTION  OF  EXPERIMENTAL  APPARATUS  AND  METHOD  OF 
ANALYSIS 


Figures  1  and  2  are  simplified  block  diagrams  showing  both  the 

field  site  recording  apparatus  and  the  laboratory  playback  facility.  A 

significant  portion  of  the  recording  apparatus  was  used  earlier  in  the  Lincoln 

wideband  noise  recording  effort  and  is  explained  in  detail  elsewhere.  The 

major  changes  are  the  following.  A  stable  crystal  oscillator  (stability  greater 
7 

than  1  part  in  10  per  day)  with  its  frequency  generator  unit  provides  a  known 
signal  level  at  a  frequency  less  than  2  Hz  away  from  the  transmitted  signal 
frequency.  This  injected  signal  level  (known  as  the  calibration  signal)  is 
used  to  calibrate  the  receiver  up  to  the  input  of  the  quadrature  correlation 
receiver  for  narrowband  linear  processing.  The  generator  also  generates 
all  the  phase  references  needed  for  the  correlation  receivers.  The  frequency 
stability  of  the  crystal  oscillator  is  checked  prior  to  each  run  by  comparison 
with  the  WWVL  standard  at  Fort  Collins.  All  of  our  processing  reported  here 
was  done  on  the  narrowband  channel  data.  [A  separate  memo  by  Evans  will 
discuss  the  nonlinear  processing  of  the  wideband  data.  ]  All  the  data  was 
recorded  on  magnetic  tape  which  were  to  be  processed  back  at  the  Laboratory. 
The  same  equipment  was  used  to  process  data  from  both  sites  at  both 
frequencies . 

The  playback  facility  at  the  Laboratory  is  shown  in  Fig.  2.  The  data 
is  played  back  with  a  speedup  factor  of  32  so  as  to  reduce  processing  time 
from  about  4  hours  to  several  minutes.  Each  of  the  signal  channels  consists 
of  a  pair  of  quadrature  correlators  followed  by  two  squarers  and  a  summer. 
The  outputs  are  then  recorded  on  chart  recorders.  The  units  marked  PAR 
amplifiers  (Princeton  Applied  Research  Model  120  Lock  In  Amplifier)  consist 
basically  of  a  multiplier  followed  by  a  single  pole  low  pass  filter.  It  can  be 
shown  that  the  equivalent  integration  time  is,  T  =  2(32)(RC).  For  most  of 
our  processing  an  RC  =  180  sec  was  used  which  corresponds  to  an  integration 
time  T  ~  3.2  hrs.  This  long  an  integration  means  that  each  4  hour  recording 
period  results  in  a  single  signal  estimate.  It  should  be  noted  that  we  have 
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Fig.  1.  Field  site  recording  equipment 


00 


Fig.  2.  Playback  facility  block  diagram  (narrowband  processing) 


identical  receiver  channels  for  estimation  of  the  transmitted  signal  and  of  the 
injected  calibration  signal.  With  the  exception  of  the  narrowband  BPF  and 
notch  filters  the  frequency  response  of  the  recording  equipment  is  flat  from 
a  few  Hz  to  about  250  Hz.  Normally  we  try  to  adjust  the  narrowband  filters 
so  that  its  frequency  response  at  the  transmitter  signal  frequency  ( f and  at 
the  calibration  signal  frequency  (f^)  is  identical.  Thus  the  received  signal 
level  (A) at  the  output  of  the  loop  antenna  is  related  to  unbiased  signal  out¬ 
puts  (denoted  by  subscript  c)  of  the  summers  and  to  the  measured  input 
calibration  signal  B  by 


A 


B 


c 


c 


(4) 


The  equivalent  field  is  related  to  the  antenna  voltage  A  by 


H 

< 0 


AJL. 
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(nD6)f 


H(P  ~  (.  007  1 7 8 )f  (5) 

where  n  is  the  number  of  turns  in  the  loop  and  D  is  the  diameter  of  the  loop. 

We  have  used  on  occasion  an  alternative  configuration  of  the  correlation 
receiver.  The  basic  difference  is  that  we  have  used  an  integrator  rather  than 
the  low  pass  filter  with  the  long  time  constant.  We  state,  based  on  experience, 
that  the  appropriate  configuration  to  use  depends  on  the  nature  of  the  noise  and 
in  general  the  results  are  the  same. 

The  important  outputs  are  also  shown  in  Fig.  2.  The  output  of  the 

2  2 

summer  is  a  biased  estimator  for  A  (or  B  )  i.  e.  , 

.  N  /T 

E[AZ]  =  A^  {1+2  -At-  j  .  (6) 

A^ 
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The  point  is  that  one  must  subtract  off  2  Nq/T  from  the  sample  value  to  get 

an  unbiased  estimate  of  A^.  Hence  this  means  that  N  must  be  determined 

o 

as  well.  The  method  of  determination  of  Nq  is  discussed  in  Appendix  A.  The 
confidence  interval  associated  with  each  measurement  A^  is  determined 
from  the  graph  in  Appendix  B.  The  confidence  interval  itself  is  to  be  inter¬ 
preted  as  a  random  variable  and  it  has  a  certain  probability  (1-c)  of  including 
2 

the  true  value  A  . 
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III.  DATA  CATALOGUE 


Due  to  a  number  of  equipment  malfunctions  at  both  the  Hawaii 
and  California  sites  a  significant  percentage  of  the  total  available  data  was 
lost.  Table  I  lists  the  numbers  of  hours  of  data  recorded  at  each  station. 

The  unfortunate  fact  is  that  the  number  of  hours  in  which  data  was  recorded 
simultaneously  at  both  sites  amounted  to  only  24  hours. 

Past  theoretical  analysis  and  experimental  measurements  provide 
estimates  of  the  expected  value  of  the  field  strengths  at  Hawaii  and  California. 
The  computed  values  for  the  direct  component  of  H^  is  shown  in  Table  2  for 
California  and  Hawaii  for  45  Hz  and  7  5  Hz  (under  daytime  and  nighttime 
conditions) . 

Tables  3  and  4  summarize  the  results  of  the  measurements  for 

California  and  Hawaii  respectively.  The  column  labelled  "H^(dbH)  -  bias 

corrected”  gives  the  estimates  of  magnetic  intensity  after  the  bias  correction 

(see  equation  6).  Noise  estimates  which  result  from  the  variance  computation 

(as  described  in  appendix  A)  are  given  for  three  time  intervals  over  the  data 

record;  at  the  beginning,  in  the  middle  and  at  the  end  of  the  data  record.  A 

number  of  noise  estimates  derived  from  the  noise  spectrum  analysis  programs 

are  shown  also.  These  spectra  were  computed  from  data  at  the  tape  beginning 

and  the  agreement  between  the  two  methods  of  noise  analysis  is  very  good 

for  the  California  data  (the  discrepancies  in  the  Hawaii  noise  evaluation 

methods  will  be  discussed  later  in  the  section).  The  last  two  columns  of 

a 

data  show  the  integrated  noise  level,  NT,  and  the  signal  to  noise  ratio  S/NT. 

A  ^2  1 

S  =  H  is  the  biased  signal  estimate  at  the  summer  output  and  Nj  the 
integrated  noise  at  the  output  of  one  of  the  quadrature  correlators.  We 
tabulate  the  biased  signal  to  noise  ratio  rather  than  the  unbiased  signal  to 
noise  ratio  because  the  first  ratio  is  the  one  to  be  used  in  determining 
confidence  intervals. 

The  factor  which  has  seriously  compromised  the  full  potentiality  of 
this  experiment  is  the  high  noise  levels  along  with  some  intermittent 
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TABLE  1 
July 


6 

7 

8 

9 

10 

Calif. 

4  hrs.  day 

4  hrs .  day 

4  hrs.  day 

4  night 

4  hrs.  day 

4  night 

4  hrs .  day 

4  night 

Hawaii 

No  Freq. 
Ref. 

No  Freq. 
Ref. 

4 hrs.  day 

4 hrs.  night 

4 hrs.  day 
4  night 

equip,  malfunc . 
4  hrs .  night 

High  Current  Nights 


13 

14  15 

16 

17 

Calif. 

2  hrs. 
day 

Frequency  Generator  Failure 

No  Data  Collected 

Hawaii 

4  day 

4  night 

4  hrs.  day  4  hrs.  day 

4  hrs.  day 

4  hrs .  day 

20 

21 

22 

23 

24 

Calif. 

1  hr.  day 

* 

1  hr.  day 

* 

data  tape  erased 
at  Laboratory 

3  hrs.  (75Hz) 
day 

tape  recorder 
failure 

Hawaii 

4  hrs.  day 

4  hrs .  day 

5!' 

4  hrs.  (75Hz ) 
day 

2  hrs  .  (75Hz)  day 

data  collection  stopped  by 
forest  fire 

Summary  of  Available  Data 
California  /Hawaii 
All  Field  Strengths  dbH 


Due  to  scheduling  error  receivers  were  set  to  receive  75  Hz  while  transmitters  were  at  45  Hz. 
As  a  result  we  have  only  wideband  data  with  no  frequency  reference. 


TABLE  2 

H  (dbH)  Estimated  Field  Strengths  (Direct  Component) 
0 


California  d  =  2.  52  Mm 

Hawaii  d  -  6.  53  Mm 

45  Hz 

75  Hz 

45  Hz 

75  Hz 

Day 

Night 

Day 

Night 

Day 

Night 

Day 

Night 

|  10  log-jj—  -  20  log  „| 

-193. 34 

-193. 34 

-193. 34 

-193. 34 

-193. 34 

-  193. 34 

-193.34 

-  193.  34 

,  „  1  =  150A 

20  log  II-  L  =  22.5Km. 

+  130.  57 

+130. 57 

+130. 57 

+130. 57 

+  130.  57 

+130. 57 

+130. 57 

+130. 57 

20  log  f 

33.06 

33.06 

37.  50 

37.  50 

33.  06 

33.06 

37.  50 

37.  50 

20  log  (cos  (o)  0  =  22° 

-.66 

-.66 

-.66 

-.66 

-.66 

-.66 

-.66 

-.66 

.  1.2  x  10"4(45  Hz) 

-10  log  crc  K7  x  i0-4(75  Hz) 

39.  21 

39.21 

37.70 

37.  70 

39.21 

39.21 

37.70 

37.70 

,  v  hn  =  75  km 

-20  log  h  h°  ,  90 

-97. 50 

-99. 08 

-97. 50 

-99. 08 

-97. 50 

-99.08 

-97. 50 

-99.08 

-  10  log  (~ — ) 
ph 

-1. 00 

-.29 

-1.07 

-.  79 

-1.00 

-.29 

-1.07 

-.79 

d/re 

10  '<*«  TUTJ^l 

.  11 

.  11 

.  11 

.  11 

.  78 

.78 

.78 

.  78 

-  1 0  log  d 

-64. 01 

-64.  01 

-64.  01 

-64.01 

-68.  15 

-68. 15 

-68.  15 

-68. 15 

~  cr(45)  =  .75  db/Mm 

'ffd  0(75)  =  * 

-1.89 

-1.89 

-3.25 

-2.  55 

-4.  90 

-4.90 

-8.  42 

-6.60 

Total  H  (dbH) 

0 

-155.  45 

-156. 32 

-153.95 

-154. 55 

-  161. 93 

-162.80 

-162.59 

-162.07 

-153.  82 

for  Ii  200A 

-160. 30 

for  1=2  00 A 

*/£ _ ) 

(V  45  Hz’ day 


=  1.26 


i£ _ i 

'Vph'  45  Hz.  night 


From  NUSC/NLL* Measurements 


75  Hz.  day 


1.28 


75  Hz. 


night 


1. 20 


0(75  Hz) 


*  Naval  Underwater  Systems  Center,  New  London  Laboratory 


day  =  1.2  9  db/Mm 
night  =  1.01 


TABLE  3 


Summary  of  California  Tapes 


Hg(dbH) 

i 

Cal.  Mon. 

Cal  Level 

Sig  Level 

H  (dbH ) 

Bias 

Noise  N  i 

Variance 

Spectral 

Integrated 

S/N  Ratio 

Date 

Level  (mv) 

PAR  (mv) 

PAR  (mv) 

tP 

Corrected 

from  °  ' 

1 

Computation 

1  Est.(No) 

Noise(LPF) 

(after  integ .  ) 

Cl 

in.* 

15.8 

1.  1 

-154.8 

-155. 3 

turn  count 

-  130.0 

Ni=N0/T 

12.  3 

-The  injected  calibration  system 

July  6 

(RC  =  80) 

(RC  =  80sec ) 

~  700 

dbH/VHl 

-167.  1 

level  B  =  (cal  mon.  level)  x 

Monday 

dbH  /^Hz  (Si 

(atten.  factor)  where  atten. 
factor  =  7.  5577  x  10*7 

(Day) 

45  Hz 

C  2 

109. 

15.  5 

.88 

-156.9 

-157. 5 

264-408 

-129. 0 

1  -129.4 

-169. 6 

12.7 

July  7 

106. 

(180) 

(180) 

694-819 

-127.7 

-168. 3 

11.4 

Tuesday 

(Day) 

1 163-1273 

-126.  1 

-166. 7 

9.8 

C3 

1  11. 

15.8 

1.  11 

-154.7 

-155.0 

325-555 

-130. 8 

•  -130.7 

-171.  4 

16.  7 

July  8 

1 1 1. 

(180) 

(180) 

830-925 

-128.  3 

-168. 9 

14.2 

Wed. 

(Day) 

1110-1280 

-127.  3 

-167. 9 

13.2 

C4 

1 11. 

16.  5 

.72 

-  158.8 

-159. 0 

140-384 

-  134. 6 

1  -135.4 

-175.  2 

16.  4 

July  8-9 

112. 

(180) 

(180) 

700-899 

-133.7 

|  -134.2 

-174.  3 

15.  5 

Wed. 

(Night) 

1145-1310 

-131.9 

-172.  5 

13.7 

C5 

115. 

16.  1 

.  87 

-156.9 

-157. 8 

270-417 

-128.  8 

-169  4 

12. 5 

July  9 

112. 

(180) 

(180) 

680-808 

-126.  4 

-167. 0 

10.  1 

Thur  s . 
(Day) 

1180-1290 

-125.  3 

-165.  9 

9.0 

C6  (200A) 

113. 

17.2 

.97 

-156. 6f 

-  156.  9f 

220-347 

-  130.  4 

UJ 

o 

o 

-171.0 

14.4 

^Current  Level  -  200A 

July  9-  10 

112. 

(180) 

(180) 

698-828 

-130. 4 

-171.  0 

14.  4 

Thurs . 
(Night) 

1130-1240 

-132. 0 

-172. 6 

15.9 

C7 

114. 

17.2 

1.  1 

-155. 4 

-156.2 

290-500 

-126. 4 

-167. 0 

11.6 

July  10 

112. 

(180) 

(180) 

715-847 

-125.  5 

-  166.  1 

10.7 

F  riday 

(Day)  

1110-1230 

-124.  3 

-164.  9 

9.  5 

C8  (200A) 

112. 

17.  0 

.  73 

-159. 0f 

-159. 8* 

135-420 

-129.  1 

1  -130.3 

-169. 7 

10.7 

*  Current  Level  -  200A 

July  10-11 

111. 

(180) 

(180) 

815-1020 

-128.8 

-  169.  4 

10.  4 

Friday 

(Night) 

1200-1368 

-129.9 

-  170.  5 

11.5 

C9 

110. 

17.  3 

.9 

-  1 57 . 4 

-158.2 

300-440 

-130. 8 

-167. 9 

10.  5 

July  13 
Monday 
(Day) 

(80) 

(80) 

650-760 

-129.0 

-  166.  1 

8.7 

TABLE  4 

Summary  of  Hawaii  Tapes 


Date 

Cal.  Mon. 
Level  (mv) 

Cal  Level 
PAR  (mv) 

Sig  Level 
PAR  (mv) 

H  (dbH ) 

<0 

■TOTOT) - 

Bias 

Corrected 

Noise  No 
j  from 

Variance 

Computation 

1  Spectral 
Est. 

Integrated 

Noise 

S/Nj(db) 

H3  July  8 
Wed. 

(Day) 

129 

64 

1. 38 

-164. 1 

-166.4 

140-386 

728-913 

980-1166 

-138. 9 
-132. 3 
-133. 8 

-136. 6 

NrNo/T 

dbH /,/Hz 
-  170. 9 

6.8 

Attenuation  Factor 
=  7.7753  x  10'7 

H4  July  8 
Wed. 

(Night) 

127 

56 

1.24 

-164.3 

-  164.9 

135-382 
730-925 
,  110-1292 

-  135. 4 

-  134. 4 
-133. 7 

-  1 37 .  5 

-176.3 

12.0 

H5  July  9 
Thur s. 

(Day) 

130 

64 

1.60 

-162.7 

-164.2 

100-365 

705-900 

[100-1280 

-  1 32 . 4 
-130. 9 
-130. 7 

-171.0 

8.  3 

H6  July  9 
Thurs. 
(Night)  

124 

55 

1.  14 

-164.7 

no  unbiased 
estimates 

140-387 

710-906 

1 195-1360 

-125.7 
-125. 5 
-  127. 4 

-128. 2f 

- 167. 0 

2.  3 

Current  Level  -  200A 
^Interference  present 

H7  July  10 

Friday 

(Day) 

129 

65 

250-390 

620-760 

1110-1219 

-127. 9 
-128.  0 
-  128. 5 

Equipment  Malfunction 

H8  July  10 

Friday 

(Night) 

126 

63 

1.23 

-  165.  1 

-167. 9* 

270-410 

678-800 

1 100- 1200 

-130. 3 
-126. 8 
-  126. 8 

-133 

-171.3 

6.2 

* 

Current  Level  -  200A 

H9  July  13 
Mon. 

(Day) 

130 

67 

1.07 

-166.6 

-169. 3 

290-410 

708-820 

1250-1370 

-133. 8 
-134. 4 
-129.0 

-134 

-  172. 9 

6.  3 

H  10  July  13 
Mon. 

(Night) 

128 

64 

.96 

-167.4 

no  unbiased 
estimates 

260-420 

690-810 

1250-1340 

-  127. 4 

-  128.  6 
-  127. 4 

-126. 5f 

-  169. 6 

2.2 

*  Interference  present 

H 1 1  July  14 
Tues. 

(Day) 

128 

60 

.71 

-169.3 

208-410 

717-840 

1184-1270 

-133. 1 
-135.5 
-133. 8 

-  173. 6 

4.3 

H 12  July  14 
Tue*. 

(Night) 

127 

58 

1. 36 

-163.8 

-165. 1 

400-700 

900-1200 

1300-1500 

-  130. 4 

-  130. 4 
-131.3 

-134 

-172.7 

8.9 

H  13  July  15 
Wed. 

(P»Yl 

132 

60 

2.05 

-159.8 

-160. 5 

300-500 

700-900 

1 100-1300 

-128. 3 
-127.2 
-126.6 

-132. 9f 

-  171.2 

11.4 

*  Interference  present 

Huge  signal  estimate 

H 14  July  16 
Thurs. 

(P»Yi 

128 

60 

1.08 

-165.8 

-167.2 

130-372 

744-935 

1 128-1297 

-129.6 
-129.  1 
-128.6 

-132. 9f 

-  174. 5 

8.7 

*  Interference  present 

HI  5  July  17 
F  nday 
(Day) 

127 

62 

.82 

-168. 4 

no  unbiased 
estimate 

300-500 

700-900 

1100-1300 

-128.  5 
-126. 5 
-126. 5 

* 

-  168. 3 

-* 

Interference 

interference  at  the  Hawaii  sites.  For  example  the  California  data  (Table  3) 
shows  noise  levels  between  -130  and  -125.3  dbH/VHz.  The  noise  spectrum 
around  45  Hz  is  flat  and  no  unexpected  interference  is  present.  These  higher 
noise  levels  have  reduced  the  S/N^  ratio  by  from  5  to  10  db  lower  that  we 
would  probably  have  observed  during  early  spring  when  the  experiment  was 
originally  scheduled.  From  Fig.  B-2  in  Appendix  B  it  is  obvious  that 
reducing  the  S/N^  ratio  from  say  20  db  to  10  db  increases  the  confidence 
interval  significantly  (i.  e.  ,  from  about  2  db  to  over  8  db  for  e/2  =  .  10). 

Table  3  (for  California)  shows  S/Nj  values  typically  between  10  and  15  db 
even  after  integration  times  of  3.2  hrs.  (the  approximate  duration  of  the  data 
record  being  4  hrs.  ).  Fortunately  the  signal  estimates,  as  poor  as  they  are, 
are  still  sufficiently  accurate  to  allow  a  number  of  positive  statements  to  be 
made  concerning  high  attenuation  rates.  This  will  be  done  in  the  next 
section. 

The  Hawaiian  data  is  even  more  seriously  compromised  by  a  combination 
of  high  noise  and  interference.  From  Table  4  we  see  the  noise  varying  from 
about  -134  to  -126  dbH///Hz.  For  some  of  the  data,  the  highest  noise  levels 
were  not  a  result  of  atmospherics  but  resulted  from  some  unidentified  time 
varying  and  intermittent  interference.  It  is  observed  that  there  are 
occasional  discrepancies  between  the  two  methods  of  estimating  noise.  There 
are  good  reasons  for  these  discrepancies.  First  the  atmospheric  spikes 
recorded  on  the  wideband  channels  were  frequently  clipped  due  to  the  tape 
recorder's  channel  dynamic  range  limitation  and  hence  the  overall  spectral 
level  is  reduced.  Clipping  however  was  not  observed  on  the  narrow  band 
channel  data  which  was  used  for  the  variance  computation.  However,  the 
variance  computation  of  noise  is  sometimes  worse  because  the  time  interval 
over  which  samples  are  taken  includes  some  intermittent  interference 
whereas  the  shorter  time  interval  over  which  the  spectral  estimate  is 
derived  may  not  include  this  intermittent  noise.  The  biased  signal  to  noise 
ratios  are  generally  between  6  and  12  db  with  four  exceptions.  The  exceptional 
cases  justify  further  comment.  Three  of  these  days  (tapes  H-6,  H- 10  and 
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H-15)  have  exceptional  high  atmospheric  noise  such  that  the  bias  term 

No  2 

(2-^-)  in  the  signal  estimator  in  effect  masks  the  signal  A  .  For  the  fourth 
case  (tape  H-ll)  the  noise  is  relatively  moderate  (~  -134  dbH/VHz)  but  the 
biased  signal  estimate  itself  is  unusually  low.  Since  the  phase  stability  of  the 
correlator  phase  references  is  known  to  be  stablejOnly  a  few  factors  can 
result  in  such  a  low  biased  signal  estimate  (provided  the  transmitter  has 
phase  stability).  There  could  be  either  a  significant  phase  shift  or 
abnormally  enhanced  attenuation  due  to  some  set  of  ionospheric  conditions. 
Whatever  the  cause  of  this  low  signal  value  there  is  not  yet  enough  additional 
data  to  support  either  of  the  above  views.  In  all  the  result  is  that  all  four 
of  these  tapes  have  no  unbiased  signal  estimate.  The  other  Hawaiian  data 
tapes  do  possess  unbiased  signal  estimates  with,  however,  unfortunately  small 
signal  to  noise  ratios. 
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IV.  INTERPRETATION  OF  THE  DATA 


From  equation  (1)  we  can  derive  a  simple  expression  relating 
the  attenuation  coefficient  only  to  the  direct  component  of  the  fields  at 
California  and  Hawaii.  Since  all  parameters  not  related  to  distance  cancel 
out  in  the  normalization  it  turns  out  that 


or 


<hJh.d 

(IVc,D 


(7) 


O' 


EW 


=  - . 870  + 


D]db 


•  [gJ iLPj 


db 


4.0 


db/Mm. 


Originally  we  had  planned  to  use  the  estimates  of  the  total  signal  strengths 
at  Hawaii  and  California  in  this  expression  and  compute  an  effective 
attenuation  coefficient  a  .  Since  we  cannot  at  present  separate  the  round-the- 
world  component  from  the  total  measured  signal  this  effective  attenuation 
coefficient  ry  will  be  in  general  different  from  the  true  attenuation  rate 
However,  for  reasonable  choices  of  and  it  can  be  shown  that  the 

expected  range  of  a  estimates  (in  the  absence  of  significant  noise)  should 
form  a  reasonable  band  of  values  around  For  example  if  one  assumes 

^EW  =  ^WE  =  *  ^  db/Mm  it  is  easy  to  show  that  a  will  fall  between  .  3 
and  1.1  db/Mm.  If  experimental  values  of  a  fall  in  this  band  it  is  reasonable 
to  assume  that  high  values  of  >2.0  db/Mm  do  not  exist.  However  as 

noise  becomes  a  significant  factor  the  expected  range  of  a  increases  to  the 
point  where  one  cannot  make  any  intelligent  judgement  of  on  the  basis 

of  the  measured  a  values.  It  will  be  apparent  that  the  extremely  poor 
S/N  ratios  encountered  in  Hawaii  prevent  us  from  determining  good 
unbiased  estimates  of  the  signal  strength.  Accordingly  the  extremely  large 
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confidence  intervals  associated  with  these  unbiased  Hawaii  estimators  are 

A 

translated  into  even  larger  confidence  intervals  on  a  .  The  point  is  that  the 
Hawaii  measurements  turned  out  to  be  sufficiently  noisy  so  that  measurements 
of  a  cannot  distinguish  between  ~  .  7  and  >2.0  db/Mm.  For  this 

reason  no  attempt  has  been  made  to  compute  an  attenuation  coefficient  from 
simultaneous  Hawaii  and  California  measurements. 

As  a  result  of  not  being  able  to  use  the  above  simple  expression  for  a 
we  are  forced  then  to  compare  experimental  signal  strengths  with  theoretically 
calculated  levels  (and  ranges)  and  make  some  quantitative  statements  on  the 
comparison.  This  puts  us  in  the  unfortunate  circumstance  of  requiring 
reasonable  knowledge  of  the  values  of  h^,  h^,  (c/V^^n  anc*  (c/V 
Some  recent  measurements  at  NUSC/NLL  have  provided  us  with  the  parameters 
we  have  used.  It  is  not  unreasonable  to  assume  at  this  point  that  the 
uncertainties  associated  with  these  parameters  may  lead  to  variations  in  the 
theoretical  signal  levels  of  the  order  of  ±1.  5  db.  Even  with  this  restriction 
in  mind  we  still  can  make  some  definite  statements  concerning  the 
attenuation  coefficient. 

It  is  evident  from  Tables  3  and  4  that  the  unbiased  signal  estimates  have 
a  considerable  spread  for  California  data  as  well  as  that  for  Hawaii.  For 
each  unbiased  signal  estimate  and  noise  estimate  a  confidence  interval  (see 
Appendix  C)  is  determined.  Likewise  a  theoretical  signal  range  for  a 
particular  ct  assumption  is  calculated  (the  range  being  determined  by  the 
magnitude  of  the  standing  wave  uncertainty  as  outlined  in  Appendix  B).  The 
next  logical  step  is  to  compare  the  theoretical  signal  range  (for  an  assumed  q?) 
to  either  the  confidence  intervals  or  to  some  grouping  of  experimental 
unbiased  signal  estimates.  We  are  going  to  consider  two  cases  of  interest 
utilizing  mostly  California  data  because  of  its  higher  S/Nj  ratio. 

Case  a  For  the  theoretical  computations  we  first  choose 

(g'd  n^EW  EW  =  «75  db/Mm  since  this  has  been  a 

commonly  used  value  in  many  studies. 
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Case  b 


For  the  theoretical  computation  we  assume  that  both 
(ffD)EW  and  are  greater  than  2  db/Mm.  For  the 

west  to  east  direction  we  assume  (q^  =  .75  db/Mm. 

Figure  3  shows  both  the  biased  and  unbiased  signal  estimates  as  well 
as  the  confidence  intervals  for  each  biased  experimental  signal  estimate. 


Table  5  lists  the  theoretical  range  of  signal  values  for  cases  a  and  b. 
The  first  row  lists  the  values  of  computed  field  strength  of  the  direct 
component  only  (values  taken  from  Table  2).  The  second  row  lists  a  weighted 
experimental  average  of  the  unbiased  signal  estimators  (H  ^)u  ^  of  the  total 
signal  where 


(H  *)  = 

O  u 


M 

£  [ 

i=  1 


<*y> 


u 


Nt 


]; 


M 

liWi 


(8) 


The  difference  between  the  theoretical  estimate  and  the  composite  experimental 
values  are  listed  in  row  3.  The  difference  under  California  daytime  conditions 
is  only  0.8  db/Mm,  suggesting  compatibility  between  experimental  and 
theoretical  results  (for  case  a).  The  other  columns  show  considerably 
larger  differences  which  will  be  discussed  later  in  this  section.  The  range 
of  theoretical  values  in  row  4  are  indicated  on  Fig.  3  as  solid  vertical  bars. 

It  can  be  easily  seen  that  the  primary  effect  of  increasing  or  decreasing  the 
east  to  west  attenuation  coefficient,  (<y^  Just  to  shift  the  solid 

vertical  bars  up  or  down  respectively.  (The  change  in  the  length  of  the  solid 
bars  being  a  secondary  effect.)  One  now  should  compare  the  unbiased  signal 
estimates  to  the  theoretical  range  of  total  signal  strength  (for  a  particular  o'). 

For  the  daytime  California  results  none  of  the  unbiased  signal  estimates 
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(dBK) 


Fig.  3.  Confidence  intervals  (California  -  Hawaii  results). 
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TABLE  5 

Theoretical  Computations 


I 


150A  f  =  45  Hz 


Case  a  <«D,N}EW,WE  =  -75db/Mm 


California 

Hawaii 

Day 

Night 

Day 

Night 

1.  H  -  Direct  Component  (see  table  2) 

40  (Theoretical  Estimate) 

-155.5 

dbH 

-156. 3 
dbH 

-161.9 

dbH 

-162.8 

dbH 

2.  (H  2)  _  Composite  Experimentally 

10  U  Measured  Value  (Unbiased) 

-156.3 

dbH 

-159. 5 

3.  Difference  =  (Theor.  Est.  -  Composite 

Exp.  Est.  ) 

+  0. 8  db 

+3 . 2  db 

4.  Range  of  calculated  signals  due  to 

standing  wave  pattern  interference 

l(HB>TlmaX&  ^Timin' 

Max.  - 155.  1 
Min.  -155.9 

-155. 9 
-156. 7 

-161.1 
-162.  8 

-162.0 
-  163.7 

Case  b  N^EW  -  db/Mm,  ^®D,  N^WE  ~  db/Mm 


California 

Hawaii 

Day 

Night 

Day 

Night 

1.  Direct  Comp.  (Theoretical  Est.  ) 

-158.6 

dbH 

-159. 5 

-170.  1 

-171. 0 

4.  Range  of  total  signal  due  to  standing 

wave  interference. 

Max.  -158.  1 
Min.  -  1  59 . 2 

-159. 0 
-160.  1 

-168. 2 
-172. 6 

-169. 1 
-  173.5 

% 


fall  below  the  lower  limit  of  range  for  (g^ew  =  2.0  db/Mm.  This  suggests 
that  during  daytime  conditions  east  to  west  attenuation  rates  are  smaller 
than  2.0  db/Mm.  Since  all  the  confidence  intervals  include  the  theoretical 
signal  range  for  (o^ew  =  .75  db/Mm  it  is  possible  that  the  attenuation 
value  is  in  fact  the  true  value,  although  the  grouping  of  the  unbiased  sample 
estimates  suggest  it  may  be  somewhat  higher.  The  likelihood  that 
(o'd)ew  <  •  75  db/Mm  also  seems  less  feasible  when  one  considers  that  the 
composite  weighted  unbiased  signal  estimate  is  about  -156.3  db/Mm. 


For  the  nighttime  California  measurements,  all  three  unbiased 
estimates  fall  considerably  lower  than  the  signal  levels  expected  for 
(o'n)ew  =  •  75  db/Mm  (in  fact  the  unbiased  estimates  group  reasonably  close 
to  the  signal  levels  expected  for  (o'^ew  =  2.  0  db/Mm).  In  addition  the 
upper  limits  of  the  confidence  intervals  do  not  reach  the  minimum  theoretical 
signal  level  for  =  .75  db/Mm.  One  must  however  be  very  careful 

not  to  jump  to  the  conclusion  that  (g^)ew  ^act  nearly  2.  0  db/Mm.  It 

must  be  realized  that  the  theoretical  signal  strengths  were  derived  on  the 
basis  of  values  of  h  and  c/V  ,  as  given  by  Table  2.  It  may  be  possible  that 
might  be  as  large  as  95km  and  ( ° jvf  as  large  as  1.2.  These  values 
would  lower  all  the  theoretical  signal  levels  (at  night)  by  about  1. 3  db.  Even 
in  this  case  a  composite  unbiased  signal  estimate  for  California  would 
indicate  an  attenuation  coefficient  closer  to  1.  5  db/Mm  rather  than  .75  db/Mm. 
The  nighttime  data  seems  then  to  indicate  an  attenuation  coefficient 
somewhere  between  .  75  db  and  2.  0  db  with  a  reasonable  probability  that  the 
value  might  be  towards  the  larger  value. 


One  can  strengthen  the  above  arguments  by  examining  the  weighted 
average  (composite)  signal  strengths  and  the  composite  confidence  intervals 
associated  with  these  averages.  The  computation  of  these  composite 
confidence  intervals  is  shown  in  Appendix  D.  Figure  4  shows  the  daytime 

The  two  values  taken  together  are  inconsistent  with  recent  NUSC/NLL  measure¬ 
ments  (which  indicate  the  values  used  in  Table  2  are  more  appropriate). 

Most  theoretical  estimates  for  (c/V  ^)j^  are  between  1.  1  and  1.2. 
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Fig.  4.  Composite  confidence  intervals. 


24 


and  nighttime  composite  unbiased  signal  estimates  and  the  associated 
confidence  limits  (for  1-e  =  .80)  along  with  the  theoretical  signal  ranges 

for  =  .75  db/Mm  and  2.0  db/Mm.  It  is  immediately  obvious  that 

under  daytime  conditions  the  confidence  limit  includes  the  range  of  values 
for  =  .75  db/Mm  and  excludes  those  for  £  2.0  db/Mm.  It  is  safe 

to  state  that  there  is  less  than  a  10$  probability  that  the  true  value  of 
is  2.  0  db/Mm  or  greater.  It  is  also  apparent  that  at  night  the  range  of 
signal  values  for  a  =  .75  db/Mm  are  well  above  the  composite  confidence 
limit.  By  making  other  less  reasonable  choices  for  h,  (c/V^)  and  cr^  it 
is  still  difficult  to  bring  the  range  down  to  the  upper  limit  of  the  confidence 
interval.  This  is  the  basis  for  believing  that  the  nighttime  attenuation  rate 
is  significantly  larger  than  .75  db/Mm. 


The  data  from  Hawaii  to  a  degree  support  (or  at  least  does  not  reject) 
the  conclusions  reached  by  use  of  the  California  data  alone.  The  confidence 
intervals  are  however  so  large  that  they  cover  both  the  signal  range  for 
a  -  .75  db/Mm  and  2.  0  db/Mm.  The  exceptions  are  the  daytime  measure¬ 
ments  on  July  13  and  14  where  the  confidence  intervals  are  below  the  signal 
range  for  o'  =  .75  db/Mm.  These  points  would  seem  to  indicate  that  the 
daytime  attenuation  rate  was  also  somewhat  greater  than  .75  db/Mm.  The 
grouping  of  the  unbiased  signal  estimates  is  so  spread  out  that  any  estimate 
of  the  probable  attenuation  rate  range  is  premature.  The  confidence  intervals 
for  both  the  day  and  nighttime  data  seem  to  indicate  that  (q^  >  .75  db/Mm. 

Since  we  have  integrated  over  such  a  long  period  (usually  about  3  to  4 
hours)  the  possibility  arises  that  a  relative  phase  shift  between  the  signal 
phase  and  the  internally  generated  reference  phase  may  occur  during  this 
extended  period.  The  result  of  this  would  be  a  signal  estimate  lower  than 
the  true  signal  level  by  an  amount  which  depends  on  the  degree  of  the  phase 
shift.  It  should  be  emphasized  that  great  effort  was  taken  to  routinely  check 
the  receivers  phase  stability.  Observations  during  the  playback  of  the  data 
tapes  revealed  no  phase  shifts  due  to  the  receiver  malfunctions.  Accordingly 
if  a  signal  phase  shift  occurred  it  must  be  due  to  changing  ionospheric 
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conditions.  Since  we  have  operated  only  when  the  entire  propagation  path 
was  entirely  in  daytime  conditions  or  in  nighttime  conditions,  one  expects 
the  phase  shift  due  to  changing  ionospheric  conditions  to  be  very  small  and 
to  effect  the  signal  estimate  to  a  negligible  degree. 

Up  to  the  present  time  we  do  not  have  sufficiently  high  antenna  currents 
to  determine  whether  such  ionospheric  phase  shifts  do  occur.  It  should  be 
mentioned  that  analysis  of  the  California  data  did  not  reveal  any  abrupt 
ionospheric  induced  phase  shift  over  that  part  of  the  path.  With  the  present 
signal  levels  one  cannot  test  whether  there  is  a  gradual  phase  shift  of  the 
signal  phase  due  to  changing  ionospheric  conditions. 
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V.  SUMMARY 


The  following  list  of  observations  is  a  summary  of  both  our  field 
experiences  and  analysis  of  the  data. 

1.  There  is  ample  evidence  that  attenuation  rates  (at  45  Hz)  are 

less  than  2.  0  db/Mm  for  propagation  over  paths  entirely  under  daytime 
conditions,  i.e.,  <2.0  db/Mm. 

2.  The  best  guess  is  that  the  daytime  attenuation  coefficient, 

^D^EW*  *s  between  .75  and  1.2  db/Mm  if  the  choices  used  for  h,  (j^  and 

c/V^  are  the  correct  ones. 

3.  At  night  there  appears  to  be  a  substantially  lower  than  expected 

signal  level.  The  best  guess  is  that  the  attenuation  coefficient  is  somewhere 
between  1.2  and  1.8  db/Mm.  The  S/N  ratios  for  these  data  (in  California) 
are  high  so  that  we  have  high  confidence  that  these  points  are  quite  accurate 
signal  estimates.  This  higher  than  expected  attenuation  rate,  is 

the  major  surprise  of  the  entire  experiment. 

4.  To  improve  upon  the  results  of  this  experiment  a  number  of 
factors  can  be  considered.  The  biggest  improvement  would  be  to  make  the 
same  set  of  measurements  at  a  less  noisy  period  (the  winter  season).  This 
should  increase  the  signal  to  noise  ratio  by  some  5  to  10  db.  A  larger 
current  level  say  300A  will  raise  the  signal  by  6  db.  It  is  doubtful  that 
integrating  over  periods  exceeding  5  hrs.  will  be  possible  since  one  wants 
to  keep  well  away  from  sunrise  and  sunset  periods  at  Wisconsin  and  Hawaii 
(about  a  2  db  improvement).  Hence  one  can  get  roughly  a  13  to  18  db  increase 
in  S/N  ratio  with  the  above.  Then  if  the  round-the-world  contribution  can 

be  separated  from  the  total  signal,  the  signal  to  noise  improvement  will 
result  in  much  better  estimates  of  a  (a  report  is  in  preparation  on  future  ELF 
propagation  measurements). 
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APPENDIX  A 


NOISE  INTENSITY  MEASUREMENTS 


Accurate  knowledge  of  the  noise  intensity  (at  the  signal  frequency)  is 
necessary  for  the  following  two  reasons. 


1.  The  estimate  of  the  signal  power  at  the  output  of  the  receiver's 

-  2  No 

summer  is  a  biased  estimator(i.  e.  ,  E[A  ]  =  A  +  2(— ^r))-  The  bias  con - 

sequently  must  be  subtracted  from  the  sample  value  A  to  give  the  unbiased 
estimate  of  A^. 

2 

2.  The  confidence  limits  on  each  estimate  A  depend  on  the  ratio 

No  2 

jj—  (where  Nj.  =  -^r).  Accordingly  a  small  confidence  interval  for  A 

requires  an  even  smaller  confidence  interval  for  Nj. 


For  purposes  of  this  experiment  we  have  set  an  acceptable  confidence 
interval  of  ±  .  2  db  on  Nj.  It  should  be  noted  that  the  best  way  to  reduce  the 
confidence  interval  in  estimating  a  signal  level  is  to  integrate  over  the 
entire  data  record  and  come  up  with  a  single  value  A^.  On  the  other  hand, 
to  reduce  the  confidence  interval  in  noise  intensity  estimation  the  scheme 
is  to  get  as  many  independent  samples  in  the  data  record  as  possible.  We 
will  defer  until  the  end  of  this  section,  the  method  which  selects  the  number 
of  samples  needed  to  provide  a  specified  confidence  interval  about  Nj. 


Figure  A- 1  shows  the  primary  method  by  which  the  noise  intensity 

NJ  was  computed.  The  output  of  a  single  correlator  (PAR)  is  sampled  at 

a  reasonable  high  rate  fsamp^e  (~  40/sec)  and  recorded  on  magnetic  tape. 

The  output  random  variable  is  denoted  as  (x  denotes  either  the  cosine  or 

sine  reference  channel).  The  sample  mean  m  and  sample  variance  a  ^ 

x  x 

are  defined  by 


—  1 


N 


m  =  E  (L  ). 
x  N  i_  i  x  x 


(a-  1) 


N-l 


N 
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i=  1 


!  'Vi  - 


—  i2 
m  J 
x J 
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Fig.  A-l.  Primary  noise  estimation  method. 
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It  can  be  shown  E[m  ]  =  A  cos  (n,  E[m  ]  =  A  sin  and 

c  s 


E[cr  ] 


No 

var  [Lx]  =  var  [Lg]  =  var  [Lc]  =  -^ 


where  T  =  2(32)RC  .  (RC  being  the  time  constant  of  the  LPF.  )  The  above 
is  true  provided  the  random  variables  (L^  are  identically  distributed  and 
statistically  independent.  [In  the  analysis  of  the  California  -  Hawaii  data 
in  general  there  was  observed  a  significant  general  increase  in  noise 
intensity  from  the  beginning  of  the  data  record  to  the  end.  Thus  we  were 
forced  to  compute  the  variance  over  shorter  intervals  when  the  total  noise 
level  was  relatively  constant  so  as  to  make  the  (L  ).  reasonably  identically 
distributed.  Values  of  the  computed  noise  intensity  at  the  beginning,  middle 
and  near  the  end  of  the  tape  are  shown  in  tables  3  &  4  in  Section  III.  ] 


The  important  point  is  that  the  sample  variance  is  computed  from  the 

N  independent  outputs  of  one  of  the  correlators.  We  use  (cr  2)  as  an 
N0  1  x  2 

estimator  for  — =- —  from  which  we  can  then  compute  the  estimate  of  N 

T  o,l 

The  details  of  these  computations  will  be  outlined  in  the  next  few  paragraphs. 

Nq  i  (as  well  as  H ^  (dbH/^/Hz  ))  will  be  interpreted  as  the  double  sided 

noise  intensity  (i.e.  ,  S(f)  =  >  -  00  <  f  <  +  «>).  Mean  and  variance 

computer  programs  are  available  for  the  Varian  620i. 


Since  the  received  signal  levels  are  given  in  terms  of  an  H  field  at 
the  input  to  the  antenna,  we  will  also  find  an  equivalent  input  noise 

intensity  field  from  the  Varian  computed  variance.  Let  us  work  backwards 

2 

from  the  Varian  computed  variance  ax  .  The  variance  at  station  4  in 
Fig.  A- I  is  then  given  by 


) 


4 


1 

(AiA2) 


(a-4) 


There  is  a  scale  factor  K  in  the  PAR  so  that  (K  being  different  for  every 
sensitivity  scale) 
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(a/),  =  “4  <°x\  • 

3  (K)4  X  4 


(a-5) 


The  next  step  is  to  relate  the  output  at  station  3  to  the  input  at  station  1.  The 
input  noise  at  station  1  is  characterized  by  white  noise  such  that 


N 


E[  nx  (fi,  t)  nL  (Q ,  t ')]  =  -|A(t-t'). 

The  output  noise  at  station  3  can  be  shown  to  be  related  by 
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For  the  tape  recorder  contribution  there  is  both  a  question  of  amplifica¬ 
tion  and  the  expansion  of  bandwidth  due  to  playback  speedup.  The 
amplification  factor  will  be  accounted  for  in  the  calibration  process.  The 
effect  of  speedup  is  to  increase  the  noise  bandwidth  by  the  speedup  factor 
F.  The  noise  intensity  in  record  time  is  related  to  noise  intensity  in  play¬ 
back  time  by  F.  The  system  gain  from  the  antenna  output  to  the  tape 
recorder  output  is  determined  from  the  calibration  process 


q  .VB"  _  rms  value  of  the  cal,  level  measured  at  the  PAR  input 
sys  B  rms  value  of  cal.  level  at  injection  point. 

Following  the  tape  recorder  we  have  been  using  an  attenuator  (with  gain 
Pj<  1).  Thus  the  noise  intensity  at  the  antenna  output  is  related  to  the  noise 
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intensity  at  the  attenuator  output  by 


N  =  N  , 

o,  ant  o,  1 


(G  P.)‘ 
sys  1 


F  . 


(a-8) 


Finally  the  noise  at  the  antenna  output  (in  volts)  is  related  to  an  equivalent 
magnetic  field  noise  intensity  by 
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Vn 


o,  ant 


N  '  ( .  007  178)f 


(a-9) 


Equations  a-4,  a-5,  a-7,  a-8  and  a-9  are  combined  to  yield 
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(a- 10) 


It  should  be  noted  that  is  also  in  rms  units  because  of  the  1 /A /2  factor 
in  the  PAR  unit.  The  results  of  this  noise  measurement  are  listed  in 
tables  3  and  4  in  Section  III  under  the  column  'noise  from  variance 
computation'. 


The  virtue  of  the  above  noise  analysis  is  that  a  number  of  noise 
estimates  for  each  tape  takes  only  several  minutes  after  the  equipment  has 
been  set  up  and  calibrated.  There  is  an  alternative  method  of  computed 
noise  intensity  which  we  have  used  on  occasion.  The  method  consists  of 
digitizing  large  segments  of  the  wideband  channel  on  playback.  The  resulting 
digital  tape  will  be  in  the  same  format  as  was  used  for  Lincoln's  wideband 
noise  study.  One  of  the  spectral  programs  produces  a  noise  spectrum  from 
which  one  can  pick  off  the  noise  intensity  level  at  45  Hz.  The  primary  dis¬ 
advantage  of  this  method  is  that  it  is  time  consuming  and  costly  in  computer 
time.  Again  on  table  3  we  have  several  spectral  estimates  made  near  the 
beginning  of  the  tape.  The  close  agreement  between  the  noise  intensity  made 
by  variance  computation  and  that  made  by  the  spectrum  estimate  is  extremely 
close.  The  advantage  of  the  spectrum  computation  is  that  it  shows  clearly 

For  the  actual  analysis  the  following  constants  were  used  A^  =  A^  =  1, 

RC=  .01  sec,  F=  32 ,  P:  = -^  and  K=  ( 10/(50xl0-3)). 
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sources  of  interference  at  frequencies  sufficiently  close  to  45  Hz  (7  5  Hz) 
which  may  affect  the  signal  estimate.  A  memo  by  Evans  [7]  will  show  a 
number  of  spectra  for  Hawaii  which  exhibit  local  interference  which  has 
degraded  some  of  the  data. 


As  a  slight  digression  we  will  outline  how  we  select  the  number  of 

A. 

samples  needed  to  satisfy  a  specified  confidence  interval  about  H^fdbH). 
Let  us  first  assume  that  the  output  of  the  correlator  (L^).  has  had  its  mean 
subtracted  off.  Denote  this  new  variable  by  For  real  time  integration 

times  of  the  order  of  .  1  sec  or  greater;  previous  wideband  noise  studies 
indicated  that  this  variable  has  a  normal  distribution.  In  particular 
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The  variance  in  this  case  is  computed  from 
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O’  itself  is  then  a  random  variable  with  a  distribution  which  can  be  shown  to 
2 

be  a  x  distribution  with  N  degrees  of  freedom  i.  e.  , 
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Hence  it  can  be  shown 
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^2 


The  mean  and  variance  of  the  random  variable  <j  are  most  easily  found  by 
determining  the  moment  generating  function  for  the  random  variable  2 . 
This  can  be  shown  to  be 
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For  independent  random  variables  we  have  (define  Z  = 
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Now  the  first  two  moments  of  the  distribution  are  found  from  equation  a- 13 
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It  is  easy  then  to  show 
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1  No  l 

From  a- 14  one  sees  that  a  is  an  estimator  of  — ^ —  and  from  a- 15  the 
variance  of  the  estimate  decreases  inversely  with  the  number  of  samples 
chosen. 
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It  is  observed  that  f  a,  (.  )  is  a  function  both  of  N  and  the  true  variance 

No  i  CT  22 

( — .  The  confidence  limits  cttt  and  aT  are  chosen  such  that 
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It  should  be  clear  that  both  q- T  and  a  are  functions  of  both  e  and  N.  For 

Lj  xj 

large  values  of  N,  Holsinger[6]  provides  some  approximate  formulas  for  the 
confidence  interval 


a  K 

=  a2  [  l  -  —  ] -2 

U  V2(N-1) 

A  K 

aL2  =  a2  [  1  +  —  e  ]  "2 
V2(N-l) 


where  K  is  defined  by 
€ 


!=— L-  r  e^dt 


'  K 


Thus  for  a  value  =  10$  and  for  a  ±.2  db  interval,  one  requires  of  the  order 
of  1, 600  samples.  With  a  sample  interval  of  25  msec  and  for  1  minute  of 
playback  data  some  2400  samples  are  normally  collected  for  the  variance 
computation. 
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APPENDIX  B  CONFIDENCE  INTERVAL  FOR  SIGNAL  ESTIMATES 

The  output  of  the  typical  receiver  structure  shown  in  Fig.  2  is  the 
2 

random  variable  A  where 


A  2  T  2  ,  T  2 
A  =  L  +  L 
c  s 


(b-1) 


The  random  variables  L  (fi,  t)  and  L  (fi,t)  are  assumed  to  be  statistically 

c  s 

independent  and  normally  distributed  with  identical  known  variances  but 
with  different  means.  The  probability  density  functions  are  therefore 


and 


1 


fL  <LC>  =  ' 

c  JZttN  j 


1 


fL  (Ls)  =  '  , _ 

s  s  yiirN^ 


(Lc~A  cos  co)^ 

2Nl 


(Lg-A  sin  (p) 


(b-2) 


2Nt 


It  is  easy  to  show  (see  Van  Trees  [7]  )  that  the  cumulative  distribution  function 
for  A2  is  related  to  the  Marcum  Q  function  by 


*2  +  a2 


P(A2  <:y2)  =  1 


y  A/Nj 


ij  IQ  (af)e 


d  l 


(b-3) 


?  ?  A  A 

P(A^  =  1  -  Q  (■ 


(b-4) 


where 


-  ys 

v  nt 


and  a  = 


_A 

N 


I 


Let  us  select  a  confidence  level  e  .  It  is  standard  practice  to  derive  the 
confidence  curves  from  the  relations 
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(b-5) 


2  =  P(a2  sV2l) 

and 

I  =  P^>H2' 

Figure  B- 1  shows  Marcum's  Q  function  with  the  curve  parameter  con¬ 

verted  to  db.  Figure  B-2  shows  the  confidence  curves  (for  fixed  levels  of  e) 
which  are  derived  from  equations  b-4  and  b-5  and  Fig.  B-l.  Using  the 
confidence  curves  is  particularly  simple.  Assume  the  confidence  level  e  is 
specified  and  that  one  measures  a  value  A^/Nj.  A  horizontal  line  is  drawn 
on  Fig.  B-2  through  this  value  of  A^/Nj.  The  two  values  of  A^/Nj  [  defined 
by  the  intersection  of  the  horizontal  line  with  the  confidence  curves  (with 
parameter  2 ) ]  form  the  limits  of  the  confidence  interval.  An  example  of 

how  to  construct  a  confidence  interval  is  also  shown  on  the  figure.  The 
confidence  interval  (A^/Nj,  itself  a  random  variable.  The  true 

value  of  A^/Nj  is  included  within  this  interval  with  probability  l-£.  From 

these  confidence  curves  it  is  easy  to  see  that  the  confidence  interval  rapidly 

^2  / 

increases  in  size  with  decreasing  values  of  A  /Nj. 

We  also  use  on  occasion  the  term  confidence  band.  To  determine  a 
confidence  band  one  first  computes  a  theoretical  signal  level  and  the  corres¬ 
ponding  signal  to  integrated  noise  ratio.  Then  a  vertical  line  is  drawn  on 
Fig.  B-2  through  the  value  of  S/Nj  intersecting  the  appropriate  pair  of 
confidence  curves  at  two  points.  The  range  of  expected  sample  S/Nj  values 
enclosed  by  two  horizontal  lines  through  the  two  intercept  points  will  be  called 
the  confidence  band.  In  brief  the  difference  between  the  confidence  band  and 
the  confidence  interval  is  the  following.  If  we  have  a  theoretically  calculated 
signal  value,  then  the  confidence  band  gives  the  range  of  sample  biased 
signal  levels  that  one  might  experimentally  measure  {(1  -  e)  100$  of  the  samples 
would  be  in  the  band)  provided  that  the  theoretical  calculated  value  were  in 
fact  the  true  value.  If  on  the  other  hand  we  are  presented  with  a  set  of 
experimentally  measured  biased  signal  levels  then  for  each  of  these  measure- 
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Pif<r) 


P(t>r)  *  Q(a,/9)  MARCUM'S  Q  FUNCTION 


1  2  3  4  5  €  7  t  9  10  H  12  43  14  19  It 


Fig.  B-l.  Cumulative  probability  distribution  (Marcum's  Q  function). 
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P(/>r)*Q 


ments  we  have  a  unique  confidence  interval.  The  statement  can  then  be  made 
that  with  probability  l-e)100$  each  interval  will  include  the  true  value  of 
signal  (or  rather  the  true  value  of  ^  /Nj). 


Fig.  B-2.  Confidence  curves  (for  Marcum's  Q  function). 
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APPENDIX  C 


ESTIMATES  OF  THE  EFFECTS  OF  THE  ROUND-THE- 
WORLD  WAVE 


It  is  now  apparent  that  the  earth -ionosphere  cavity  can  have  a  standing 
wave  field  pattern  with  significant  nulls  and  peaks  provided  that  the 
attenuation  rate  (in  at  least  one  direction)  is  small.  This  is  normally  the 
case  for  propagation  in  the  lower  ELF  domain.  Under  the  conditions  of  a 
world-wide  isotropic  and  homogeneous  ionosphere  of  constant  height,  the 
total  field  may  be  decomposed  into  a  direct  component  and  an  around-the- 
world  component  given  by 


and 


<h„>d 


=  e 


-O' 'd 


jk 


o  V 


ph 


(c-1) 


(H0lR  =  e'  ff'(2Tre’d)  e  Jk°  V  (27rre'd)  e  j  v/2 

FR 


(c-2) 


where 


F^  =  F 


ILf 


D  R  2?]h 


COS  (f) 


/ 


2itu 


V4 


r  C  O  Lr 
e  e  V 


ph 


-)sin  0 


(c-3) 


and 


j  277  r-d 

.  a  ,  e  v 

sin  0  =  sm —  =  sm  ( - )  . 

r  r 

e  e 


If  the  ionosphere  is  not  isotropic  and  homogeneous,  then  a  common  practice 
has  been  to  replace  the  exponential  dependence  by 


d 

a 'd  -  f  o-D(r)dr 
o 


I 


o 


dr 


(2TTre-d) 

a'(2 7rre-d)  -  J  (yR(r)dr  ;  etc. 


40 


resulting  in  equations  2  and  3  in  Section  I.  It  must  be  mentioned  that  at 
present  there  is  no  theoretical  justification  for  the  simple  replacement  of 
exponential  terms.  However  it  tends  to  accentuate  the  nulls  and  peaks  of 
the  standing  wave  pattern  and  hence  provides  good  bounds  for  the  true  value 
of  total  signal  strength  (H^)rp.  Figure  C- 1  shows  curves  of  [(*^-0  |  and 

|iHoj.R  j  for  various  constant  values  of  and  a In  absence  of  phase 
information  one  can  compute  the  maximum  and  minimum  values  of  from 


I<HJtUx  KHJpl*  I'HJr 
"  I(h«>d 


(c-4) 


HHm>T  Imin 
Figure  02  is  a  graph  of 


|  (H  ) 


(n  T  1  min 


I'n.'o 


l<HJp 


and 


PVd 


Imax 


i«h«!r 


I<Vd 


We  use  these  curves  in  the  following  ways: 


versus 


iiiysd 

l(H„»D  I 


c-  5) 


(1)  If  the  attenuation  coefficients  and  are  constants  one  can 

use  curves  similar  to  those  in  Fig.  C- 1  to  determine  |(H  )^|/F  and 

|(H  I/F  as  well  as  the  ratio  l  .  If  these  are  known  then  the 

<n  R  I^dI 

maximum  and  minimum  values  of  (H  )_  can  be  found  from  Fig.  C-2. 

<P  L 

(2)  If  the  attenuation  coefficients  and  vary  with  direction  and 

position  then  one  must  use  the  full  equations  2  and  3  of  Section  I.  Then  with 

the  use  of  Fig.  C-2  one  can  compute  the  maximum  and  minimum  values  of 

(H  ) 

<0  T 
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lap)  |d/aa(*H)| 


Fig.  C-l.  Normalized  direct  and  round-the-world  field  strengths. 
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Fig.  C-2.  Maximum  and  minimum  standing  wave  interference  effects . 
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APPENDIX  D 


DERIVATION  OF  A  COMPOSITE  CONFIDENCE  LIMIT 
FOR  A7 


In  the  text  we  have  defined  the  unbiased  composite  field  strength  as 


M  (A.  ) 

T.  1  U 


.2  A  i=l 
A  = 


N. 


2  A  ? 
(A^  =  H  *  ) 
O 


(d-1) 


M 


1 


?  .NT 

1=1  1 


where  INh  is  the  ith  integrated  noise  estimate  and  (A^  is  the  corresponding 
unbiased  estimate  for  A^.  It  is  noted  that  the  weighting  gives  greater  weight 
to  those  samples  which  exhibit  the  higher  S/N  ratios.  The  1ST  can  be 
determined  to  a  very  high  degree  of  accuracy  and  for  purposes  of  computing 
a  composite  confidence  limit  these  values  will  be  considered  known  constants. 
The  mean  and  variance  of  A^  then  are 

M 


E[A2]  =  1-1  1 


M  i  -> 

A2 

1=1  1 


M 


I 


M 


=  A‘ 


1=1  1 


1 


.E.± 
1=  1  1 


and 


Wjz™  l<AXl 

var  (A2)  =  -1=1  1 


(I  — )2 
'  N.  1 

l 


(d-2) 


For  the  receiver  structure  exhibited  in  Fig.  2  it  can  be  shown  that 


A2  AZ  2  2 

var[(Apu]  =  var  [A/]  =  4  (N.*  +  N.A*) 


(d-3) 


for  A  /N^  >10  we  can  make  the  approximation 


var  [  (A.2)  1  =“  4A‘'N.  and  hence 
1  1  uJ  1 


(d-4) 
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N 


4 A 


var  (A  )  = 


.  .  N. 

i=l  1 

N  ,  2 

(  L  — ) 

i=l  l 


4 A 


N  J_\ 
(  r  n. 

ui  1 


2  2 
=  KTA^ 


(d-  5) 


where  K  = 


M  _1_ 
L  N. 


i=  1 


The  random  variable  A  can  be  assumed  to  be  Gaussian  distributed  in  the 

2 


vicinity  of  the  mean  A 


f(A2;  A2)  = 


2  2  2 
(a*  -  A*r 

2(AK)2 


yiir  ak 

The  confidence  curves  are  derived  from  the  relationship 


(A2)u=y+A2 


f(A2;  A2)dA2  =  1  -  e  . 


(a2)l=-y+a2 


Let  us  first  provide  a  normalization  by  setting  y  = 
reduces  to 


AK  -  Y 

J  t— IT -  e  2  ]  dy  =  1  -  e  . 

.  X-  'J2* 

AK 


A2  -  A2 

AK 


(d-6) 


(d-7) 


Equation  d-7 


The  upper  and  lower  bounds  are  both  determined  from  the  following 

t— = 


0(Sk>  =  *  f 


±  CD 


y_  /2tt 
AK 


2 

_ 

■  2  J  *  I 


(d-8) 
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Fig.  D-l.  Composite  confidence  curves. 
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e 

2 

.  10 

.  05 

.  02  5 

Y  -  a 

AK  P  C 

B(  =  1-28 

B  =  1.64 

0f  =  1-96 

Thus  the  upper  and  lower  confidence  curves  are  derived  from 


(AZ)U  =  A2  +  (AK)  (d-9) 

(A2)l=  A2  -  0e(AK)  . 

Figure  D-  1  shows  typical  composite  confidence  curves  for  the  averaged 
daytime  and  nighttime  unbiased  signal  estimate.  The  confidence  intervals 
(at  c/2  =  .  1)  for  the  composite  values  (A^)^  and  (A^)^.  are  also  shown. 
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